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Abstract 

Quantization and asymptotic behaviour of a variant of discrete random walk on integers are 
investigated. This variant, the Syk walk, has the novel feature that it uses many identical quantum 
coins keeping at the same time characteristic quantum features like the quadratically faster than 
the classical spreading rate, and unexpected distribution cutoffs. A weak limit of the position 
probability distribution (pd) is obtained, and universal properties of this arch sine asymptotic 
distribution function are examined. Questions of driving the walk are investigated by means 
of a quantum optical interaction model that reveals robustness of quantum features of walker's 
asymptotic pd, against stimulated and spontaneous quantum noise on the coin system. 
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I. THE eyk QUANTUM RANDOM WALK ON Z 

Recent research activity on the topic of the so called quantum random walks has achieved 



a number of interesting results and extensions of t 



le usual notion o: 



the following articles and reviews 
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random walk (see e.g 
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13|). Our present 



study investigates further quantum walks and their asymptotic statistical behaviour. Let us 
consider a walker system with Hilbert space Hyj = span{\n > \ n E Z}, and a coin system 
with space He = span{\+) , |— )}. Let the evolution operator be of the form 

y = (p+ ® ^+ + p_ ® E^)u ® 1. (1) 

This acts on the tensor product space He ® Hw Here [/ is a unitary operator act- 
ing on the coin Hilbert space and P^,P- are the projection operators in the distin- 
guished basis {|±)} of He- Also E± are the right/left one step operators in the distin- 
guished basis {\n) \ n E Z} of Hw As they are commuting they share the eigenbasis 
< 10 >= ^ e"*"*^!?! > 10 G [0, 27r] > , of orthogonal elements viz. < 0|0' >= ^5(0 - 0'); 
which satisfy the eigenvalue equations -£"±10 >= e^*'^|0 > . 

The evolution operator V acts on the state |0) to give V |0) = V{<f)) |0), where V{(f)) = 
(e*'^P+ + e~*'^P_) U = e^'^'^^U. Here is a unitary operator that acts on H^. Suppose 

now that we consider the Syk model in which the one-step evolution of the walker 

density matrix p° is given by = eyk{p^) = Trn^iV^pc ® p°V^^^), where pc is the coin 
density matrix. Here Syk is a completely positive trace preserving (CPTP) map that 
acts on Huj H^, however, since it has the property that it maps density matrices to 
density matrices we can reduce its action to the convex subset of density matrices D{Hyj) C 
Hu, ® H^. Suppose that p'^ = J J p(0, 0')|0 >< 0'|- We have that eyk{\(f) >< 0'|) = 



Trn^ {V^H<P')V'^{<l>)Pc) 10 >< 01 = ^(0,0010 >< <P'l hence p' = //p(0, 00^(0, 0O|0 X 



0'|. The n— step evolved walker density matrix is p"" = j j p(0) 0')"I0 01) 



where v4(0, 0') is the characteristic function of the quantum random walk. 

The position observable L\m >= m\m >, m E Z, and its positive powers lend themselves 
to study the statistical moments of the quantum walker after n steps. To emphasize the 
connection with the classical walk we also introduce a sequence of classical random variables 
L^"'' over the common sample space Z, that correspond to the positions of the classical walker 
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after n steps, with probability distribution P{Lr' — m) — Tr(|m)(m|p") = (m|p"|m). Then 
we obtain for the statistical moments 

(2) 

A further study of the asymptotic behavior of the first and second moment leads to the 
following 

Lemma 1 The mean position of the quantum random walk is of the form /in — {L)n — 
Kin + Hq , while the variance is of the form — {L^)n — (L)^ — K2n^ + K2,n + a^, where 
Ki,K2,K3 depend on the initial coin density matrix, the initial walker density matrix and 
the tracing scheme, however they are independent of the number of evolution steps taken, and 
fiQ and o"q are the mean and variance of the initial position distribution. In particular, if 
the quantum walker is initially in the at the position 0, that is po — \0 >< 0\, then /in — Kin 
and — K2n^ + K^n. 

Remark on L'^— quantization: Let us consider the general QRW as introduced above, in 
the particular case where U — 1, we have that £yfe(p°) = pc++E^p^ E'^^ + pc- -E'tp^E'^]^. 
This admits the interpretation that the walker shifts its position k steps to the right (left) 

with probability Pc++ (Pc )■ So if we are given a classical random walk on Z with nearest 

neighbor transition probabilities p, 1 — p, we can select a coin density matrix whose diagonal 
elements are these transition probabilities in the preferred coin basis, and view in this way 
the classical random walk as an Syk quantum random walk with U = 1, and step size k. A 
natural way then to quantize this classical walk is a unitary choice so that U ^ 1. This 
amounts to a continuous deformation of classical walk that may result into a new walk with 
novel quantum features, as its known. Also since U acts only on coin space, and we also need 
to insert a consistent coin density matrix, this U — quantization procedure as may be called, 
should be understood as a quantization of the coin of a classical random walk. Generally 
the so obtained [/— quantized random walks may admit solutions that employ off diagonal 
walker density matrices, and these are considered to be genuine quantum mechanical random 
walks. In fact the U— quantization is not restricted to walks on integers, but it can be used 
more generally to quantized classical walks on other fields of numbers and general lattice 
systems. 
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II. ASYMPTOTIC BEHAVIOUR OF QUANTUM RANDOM WALK 



Next we study the behaviour of the walker when the number of steps n, is large. In this 
case we have that 

j d(t>p{<t^,(t>)[d4>A{(^A%=^ + 0{n^-') = ^J dM<l>,<l>)K<l>r + 0(n^-'). (3) 



Here we have defined the asymptotic characteristic function h{(j)) of the walk as 

h{<i>) = -I [d^A{<i>,<i>%^^ = ImTrHXV'Kct^) {V\<i>))' p,). (4) 

For the Syk model it reads 

h{<t>) = Tth^ [{a + V\(^)aV{ct>) + • • • + V^^-\(^)aV^-\c^))p,] (5) 

where a — a^U is, & rotated Pauli matrix. 

As the following theorem states the sequence of probability measures {^Tr{\m) {m\p^) — 

r(") 

P{-^ — n — 1,2, ...} has all its moments converging to the moments of y = h{(f), as 
n goes to infinity, where (j) is stands for a random variable (rv) on the circle with measure 

r(") 

Theorem 2 The sequence of classical random variables -^,n' — 1,2, ... , corresponding to 

the sequence of quantum observables, converges weakly to the random variable Y = h{(p), 
where (f) is now a random variable with values on the circle and measure ^p((f),(/)). The 
probability distribution function for Y is given by the formula 

yi<h{(t>)<y2 * 2/1 

where it is assumed that locally in the interval [yi,y2] the function h admits a number of 
local inverses, labelled by i. 

Remarks: i) The dual e*{X) of a CPTP e{p) = J2i ^ipA] operating on a observable X and 
a density operator p respectively is obtained to be £*{X) = '^■AlXAi ii) The asymptotic 
characteristic function has a finite Fourier expansion, hence it is locally invertible except from 
a set of measure zero, if the function is not a constant. Further the inverse is differentiable 
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except from a set of measure zero. Hence the asymptotic distribution of Y given above 
holds as long as h is not a constant. A somewhat similar theorem is obtained in ^| for 
the original QRW with two essential differences. First, in Q] the sample space of rv (j) 
was enlarge by taking a product of the circle with a set of two points. Second, the initial 
walker-coin state information was encoded in the measure while in the present case only the 
initial walker state information is encoded in the corresponding measure. 

It is now reasonable to ask what happens if h is constant function. This is answered by 
the following lemma: 

Lemma 3 // the asymptotic characteristic function h is constant then cr,^ = {L'^)n — {L)^ = 
K^n + CTg, which means the quantum random walk spreads classically. 

This can be proven simply by considering the asymptotic expansion of the moments. 
As an example, let |+ > be the initial coin state, and choose U = e**'^^. The walker state 
is initially |0 >, while the model is taken to be the Sy^- In this case it turns out that 
h{4>) = — cos(20), while p{4>,(p) = 1. 

Remarks: i) The range of values of h{(p) is the interval [—1, 1]. Given that the possible 
values of walker position after n steps extends on the interval [—2n,2n], the normalized 
position asymptotic range of values is expected to be the interval [—2,2]. This means that 
the asymptotic position distribution experiences a sudden cutoff at half the spread one would 
expect. This cutoff has also been observed with different quantum evolution schemes 
13, 0, ^0 Since takes values in the interval [0, 27i], there are four relevant inverses 

of h with domain the interval [—1,1], and range in [0,27r]. All these inverses satisfy the 



relation h'(h^ (y)) = ±2a/1 — y^. The resulting asymptotic distribution is 

P{yi <Y<y,) = - [ -7^=dy. (7) 

yi 

It is worth mentioning that this is the same distribution as the one obtained in in a 
different context. 

In fact a more general result about the asymptotic distribution can be proved along the 
same lines as the following theorem states: 

Theorem 4 Whatever the initial coin density matrix and whatever the unitary coin reshuf- 
fling matrix for the ey2 quantum random walk, if the asymptotic characteristic function h{(p) 



is not constant then the normalized random variable = is distributed according to 
the distribution of eq. (0). 

This theorem is to be understood as a quantum ey2 version of the De Moivre-Laplace 
theorem for the normal approximation of the binomial distribution. In fact the correspond- 
ing theorem one would obtain, if the quantum system followed the ey evolution, would be 
precisely the De Moivre-Laplace theorem. In the case of the 5yn+i QRW the asymptotic 
characteristic function is of the form = /i — J2m=i cos2m(0 + am), where and 
am depend on the coin state and the reshuffling matrix U. If we define Ym = cos2m(0 + am) 
then it is not difflcult to show that it is distributed according to eq. ((Tj) for any m. Hence 
we obtain the following proposition: 

Proposition 5 // the asymptotic characteristic function is not constant then the normal- 
ized asymptotic position random variable Y is written as Y = + Ylm=i ^mXm, where each 
Ym is distributed according to eq. with ay = |X]m=i^m- (^^ should be noted however 
that the rv's Ym are strongly correlated) . 



III. CAVITY DRIVEN QUANTUM RANDOM WALK 

To probe the effect of initial coin state upon the long time behaviour of quantum walker 
system we introduce an operational way to exercise control on that coin state. If the coin, 
which is taken to be a two level atom, is prepared through its interaction with a single 
electromagnetic (EM) mode optical cavity on resonance and if the coin-cavity interaction is 
assumed to be of Jaynes-Cummings model (JCM), (see and references therein, also for 
the case of more general models), then its dynamic state is described as 

euipc) = TrjU{t)ipc (g) Pf)U{ty = SoPcS^ + S^pcSl (8) 

The CPTP map Su, describes a stimulate transition of the coin quantum system due to its 
interaction with the cavity mode, and is determined by the unitary evolution U{t) of JCM, 
and the initial field state which is taken to be a pure state pf = |/)(/|, where / stands 
for the cavity photon number. The state of atomic coin system after crossing the cavity is 
described by the reduced density matrix given above, where t stands for the crossing time. 
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The general case of an initial sharp number state |/) = |r), r = 0,1,2,..., leads to the 
reduced coin density matrix 



eu{pc) = AipcA\ + A2pcA + ^Pc4 



(9) 



with generators 



^1 



A, 



cos(AtVrTT) 













sv[l{Xt^/r) 




cos(Atv^) 



sin(At\/r + 1) 



(10) 



The trace preservation of this map requires that A\Ai + A|A2 + A\A^ — 1. 

Consider the special case pc — |c) (c| where |c) = cosx|+) + ^sinx|— ), that leads to 
symmetric walk about the origin. After the coin crosses the JCM cavity, its state becomes 

£u{pc), where 



£[/(Pc) ^ ^1 + ^ sin(2x) cos Xy/r + l)tj cos ( X^/rt) a-i 

cos -Jri) sin^ x 



1 

+2 



0-3 • 



(11) 



Suppose now the prepared atomic coin is used to drive a quantum random walk according 
to the £yi model for U — e*^'^^. In this case the asymptotic characteristic function reads 



X 



— COS cos^ X + cos (2Xty/r) sin 

sin(2x) cos ^ XWr + 1)J cos ( Xt^/r) sin(20) 



cos (20) 



(12) 



To evaluate the limit probability distribution function, we rewrite the last equation as 
h{(j);Xi't) — C{x-,t)cos[2(f) — A(x, f)], where we have introduced the functions A{x-,t) — 
- cos {2XtVr + 1) cos^ x + cos (2At-v^) sin^ x, B{x, t) = sin(2x) cos ( XtVr + l) cos ( Xt^/r) , 
from which we define the two functions C(x, t) = \/ A{x, t)"^ + B{x, ty and tan A(t) = . 
If we define the rv y = /i(0; t), then the normalized asymptotic walker position distribution 
becomes 



P{y;x,t) 



n^C{x,ty-y''^ 



-l<y<l. 



(13) 
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This pdf depends on the crossing time t and on the initial coin state through x- The mean 
and standard deviation derived from the above pdf are respectively 

^ = lim(-)„= / h{cP;x,t)^ = 0, (14) 



a{x,tr = lim(f^) )„ = / hic^;x,tr^ = ^^. (15) 







It is useful to observe at this stage that the distribution of eq. (fT!^ . is robust 
under changes in unless it happens that C(x,t) = 0. In this case the distribution 
collapses. Inspection of the functions A{x,t), B{x,t) and C{x,t) given above reveals that 
if X = {0, f ,7r, namely if |c) = -|+), -«|-)} respectively, and t = {§^=, 

^ 4V^"" ' 4At/r^ ' ^ 4A^"" } i k E Z, then C(x,t) = 0. The specific relations among the four 
coin states and interaction times as given above result into only two different pairs of coin 
density matrices and interaction times namely, Pc = \+){+\, t = Iav^TT ' P'^ ~ 
t = for which C(x,t) = 0. If either of these two conditions occur we say that a 

resonance condition takes place between the field and the two level atom. In this case, 
< >n ^ n, so we loose the quadratic diffusion time speed up, characterizing the quantum 
random walk. In such a case the asymptotic behaviour of the standard deviation agrees 
with that of a classical random walk. More precisely what happens is that in all the above 
cases, the exited coin is in the maximally mixed state eu{pc) = ||+)(+| + = If 

this coin system is used in V"^ QRW, then the final one-step density matrix for the walker 
system becomes 

11 1 

Pw £v^{pw) = ^Pw + -EIpwE^_^ + -E^pvy^- • (16) 

Due to the last equation, if pw is diagonal initially then so is finally. Hence we really have 
a classical one-step transition that leads to Gaussian statistics for large n, once we normalize 
L to . That implies that on resonance the walk becomes fully classical. This analysis 
makes obvious the fact that a judicious choice of the initial coin state permits us to tune 
the interaction time in a Jaynes-Cummings cavity, so that we have a quantum to classical 
transition on the asymptotic behavior of the ey^ quantum random walk. This conclusion 
makes the quantum optical experimental investigation of this idea worthwhile. 
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IV. SPONTANEOUS EMISSION IN COIN SYSTEM 



What we are going to study in this section is the effect of a spontaneously emitting coin 
system on the QRW evolution. We will assume that the coins that come into contact with 
the walker system are corrupted by spontaneous emission from the state |+) to the state |— ). 
The state |— ) is assumed to be decay stable, and state |+) is metastable and has probability 
7 of decaying to the state |— ). Let the initial coin state be Pc = |— ) (— | + pi |+) (+|, 
then the effect of spontaneous en.t.ion is to n^odtfy tMs denshy n^atnx to p' ^ (p„ + 
IPi) \~) + (1 l)Pi 1+) (+1- This effect can be captured by the generators |17|. 

r„|->H-> T.H^o ^^^^ 



In terms of these generators the Sse CPTP map that transforms to is 

Pc = esEip'c) = ToP%Tl + T,p%Tl (18) 

Of course TqTq + tITi = 1. In fact such a spontaneous de-excitation can happen as result 
of a stimulated transition as e.g in the case where the coin passes through an initially empty 
quantum optical cavity, (see the case of JCM previously), in which case 7 = sin^ At, and 
Tq = 5*0 (t), Ti = 5*1 (t), where t stands for the time spent by the coin inside the empty cavity. 
Indeed for the particular initial vacuum state |/) = |0), we obtain the map 

Pc = ^uip'c) = esEipl) = S.pcSl + S,pcS\ (19) 
where the so called Kraus generators of that map are 

So{t)= I " I , S,{t)= I ° ° I , (20) 




sin(At) 

and satisfy the property S^Sq + S\Si = 1. Since in this case no EM field is in the cavity, 
Eu corresponds to spontaneous decay of the excited atom state (as opposed to stimulated 
decay if |/) is not |0)), through the emission of a photon. 

Suppose now that we have Sy^ evolution with U = e*?'^^, and that the coin state is not 
a classical mixture of states but rather that the coin starts up in the state |c) = cosx|+) + 
ismx\—), Pc = |c) (c| , then eq. (fTT|) tells us that 

Pc = ^1 + ^yi^sin 2x^2 + ^[(1-27) cos^x-sin^xjas. (21) 



Since in this case A(x, 7) = —[(1 — 27) cos^ x ~ sin^ x] ^-^id B{x, l) — Vl ~ 7 sin 2^, we have 
that C(x,7) = ^/[(l — 7) cos 2% — 7]^ + (1 — 7) sin^ 2%, and the asymptotic characteristic 
function is ^(0; x, 7) = A(x, 7) cos 20 + 7) sin 2(f). This gives us that the mean position 
is zero and that the standard deviation is cr(x, 7) = '^'^^^ ■ Note that even in the case of very 
strong decay (7 = 1) the asymptotic mean position of the walker remains zero. This is due 
to the symmetrizing effect of U. However the standard deviation of the position depends 
exphcitly on the decay rate. Nevertheless it is not monotonic in the decay rate. Since 
(T^(x, 7) = |[1 — 47(1 — 7) cos^ x], when 7 = or 1, we have = |, and this is maximum 
spread. Furthermore the spread decreases until 7 = | and then increases again. So we can 
say that the effect of spontaneous emission is to decrease the spread of the walker position 
distribution, however the minimum spread is achieved when the probability of decay is \. 



V. CONCLUSIONS 



Here we have analysed the asymptotic behaviour of the Syk quantum random walk on 
Z. Firstly, we see that generically the standard deviation of the position distribution in- 
creases linearly with the number of steps, which is a quadratic speedup over the classical 
random walk. However there are resonant quantum random walks in which the asymptotic 
characteristic function is constant, which spread classically. Off resonance the normal- 

r(") 

ized position random variables are shown to converge weakly to the random variable 
Y — h{(f)), where (p, is a random variable whose distribution depends on the initial walker 
density matrix. Furthermore, the asymptotic distribution of the ey2 model, off resonance 
is obtained and it is found that if the position random variable is normalized, then the 
asymptotic distribution assumes a universal form given in theorem 4. Secondly, it is shown 
how to drive the asymptotic behaviour of the ey2 quantum random walk through the use 
of a Jaynes-Cummings model interaction that prepares the coin system. It is found that 
by a judicious choice of the initial coin state, it is possible to drive the quantum asymptotic 
behavior of the walk to classicality by tuning the time spent by coin to cross the Jaynes- 
Cummings cavity. Finally, the effect of spontaneous emission of the coin is considered and 
it is found that in the particular case of £y2 evolution with U = e*^'^^, it preserves the 
symmetry of the walk, and it actually decreases the spread in the walker position. 
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